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The  structures  of  Ca6oMgi5Cu25,  Ca6oMg2oCu2o  and  CasoMg25Cui5  metallic  glasses  have  been  investigated 
by  neutron  and  X  ray  diffraction.  The  correlation  functions  show  a  peak  manifold  in  the  region  2.2  4.5  A, 
arising  from  the  various  atom  pairs  for  these  glasses.  The  results  show  clearly  that  there  are  contacts 
between  the  solute  atoms  (Cu  and  Mg)  which  agree  with  a  simple  estimate  of  the  maximum  solute  atom 
fraction  beyond  which  solute  solute  contact  becomes  topologically  necessary.  The  Cu  Mg  and  Cu  Ca 
distances  are  consistent  with  the  sum  of  covalent  radii,  whereas  all  other  interatomic  distances  are 
consistent  with  the  sum  of  metallic  radii.  The  neutron  and  X  ray  diffraction  correlation  functions  were 
simultaneously  fitted  with  a  series  of  symmetric  peaks  to  obtain  coordination  numbers  and  interatomic 
distances,  but  only  the  results  from  the  first  two  peaks,  Cu  Cu  and  Cu  Mg,  are  of  reasonable  reliability. 
The  Percus  Yevick  approximation  for  binary  hard  sphere  systems  has  been  used  to  simulate  the  results. 
This  shows  that  the  results  of  fits  to  the  Cu  Cu  and  Cu  Mg  correlations  give  a  reasonable  description  for 
atoms  which  are  in  close  contact,  but  do  not  include  an  additional  contribution  at  longer  distance  due  to 
atoms  which  are  almost  in  contact.  The  fits  to  the  Cu  Ca  peak  are  strongly  affected  by  overlaps  with  both 
the  longer  distance  peaks  used  in  the  fit,  and  the  broad  trailing  edge  in  the  distributions  of  interatomic 
distances  that  is  not  taken  into  account  by  fitting  with  symmetric  peaks.  Final  results  show  that 
copper  has  a  total  of  6.5  neighbours  which  are  in  close  contact,  but  has  a  total  coordination  number  of 
about  12  13,  when  atoms  which  are  almost  in  contact  are  included. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Metallic  glasses  have  a  wide  range  of  potential  applications  due 
to  their  magnetic  properties,  high  yield  strength,  low  elastic 
modulus  and  corrosion  resistance.  The  earliest  reports  of  metallic 
glasses  were  for  compositions  which  required  extreme  cooling 
rates  to  form  minute  quantities  [1],  but  since  that  time  many 
compositions  have  been  reported  which  can  be  quenched  at 
a  slower  cooling  rate  to  form  bulk  amorphous  materials.  This  has 
allowed  the  unique  properties  of  amorphous  metals  to  be  utilised 
in  the  development  of  a  wide  range  of  products,  including  small 
scale  devices,  such  as  micro  electro  mechanical  systems  (MEMS) 
[2,3].  Despite  a  wide  range  of  current  and  potential  technological 
applications,  a  detailed  understanding  of  the  link  between  the 
structures  of  glasses,  the  glass  forming  abilities  and  the  physical 
properties  they  exhibit  has  not  been  reached  and,  although  models 
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for  the  structure  of  binary  glasses  have  been  proposed,  the  struc 
ture  of  more  complicated  amorphous  systems  is  still  not  well 
understood. 

A  number  of  structural  studies  of  bulk  metallic  glasses  have 
been  carried  out  using  both  X  ray  and  neutron  diffraction,  as  well 
as  by  applying  modelling  methods  such  as  Reverse  Monte  Carlo 
(RMC)  and  molecular  dynamics  (MD).  Significant  progress  has 
been  made  in  understanding  the  structure  of  binary  metallic 
glasses  using  these  methods,  and  the  results  indicate  that  the 
packing  of  atoms  in  these  materials  is  not  random,  but  is  strongly 
influenced  by  chemical  interactions.  In  the  case  of  NigiBig,  for 
example,  there  is  a  strong  preference  for  unlike  atom  bonds,  and 
there  are  no  B— B  nearest  neighbours  at  all  [4].  These  studies  have 
resulted  in  the  formulation  of  the  efficient  cluster  packing  (ECP) 
model  [5].  In  this  model,  the  simplest  case  is  one  for  which  the 
solute  atom  (B  in  the  case  of  Ni8iBig)  is  relatively  dilute  in  the  alloy, 
and  all  the  nearest  neighbours  are  solvent  atoms  (Ni  atoms).  In  the 
literature,  there  are  many  examples  of  glasses  for  which  this 
solute  lean  model  is  applicable,  including  NiglBig  [4],  F8oB2o  [6], 
Ti84Sii6  [7]  and  C080P20  [8].  However,  many  of  the  bulk  metallic 
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glasses  which  exhibit  the  modest  cooling  rates  and  large  critical 
diameters  necessary  for  bulk  applications  are  complex  systems 
containing  three  or  more  atom  types,  and  with  solvent  concen 
trations  that  are  too  low  to  allow  each  solute  atom  to  be  sur 
rounded  only  by  the  solvent.  Instead,  following  the  idea  of  the  ECP 
model  for  solute  rich  glasses  [9],  solute  atoms  must  substitute  onto 
solvent  atom  sites.  It  is  thought  that  the  more  complex  structure 
this  produces  is  the  cause  of  the  increased  glass  forming  ability 
(GFA)  of  these  materials  because,  coupled  with  the  high  density, 
and  therefore  low  free  volume  of  a  bulk  metallic  glass  [10,11  ],  it  is 
much  harder  for  the  atoms  to  rearrange  to  form  a  crystalline 
structure  [3].  These  solvent  anti  site  defects  can  also  increase  the 
number  of  unlike  atom  bonds,  further  increasing  the  stability  [9], 
Senkov  et  al.  [12—16]  have  published  several  papers  on  the  GFA 
and  physical  properties  of  the  glass  systems  CaxMgyCu(ioo-(x+y)), 
CaxMgyZn(ioo-(x+y))  and  CaxMgj,ZnzCu(10o-(x+y+z)).  These  studies 
demonstrate  that  the  GFA  of  the  system  is  best  for  compositions 
where  the  solvent  (Ca)  accounts  for  approximately  half  the  atoms. 
For  example,  the  composition  with  the  best  GFA  for  the  ternary 
system  containing  copper  is  Ca5oMg22.5Cu27.5,  where  the  critical 
thickness  is  10  mm.  This  reduces  rapidly  as  the  amount  of  solvent 
(Ca)  is  increased,  with  samples  containing  >60  at.%  Ca  having 
a  critical  thickness  of  no  more  than  4  mm,  reducing  to  <1  mm  for 
two  compositions  with  70  at.%  Ca.  Furthermore,  the  GFA  is 
improved  when  the  number  of  Cu  and  Mg  atoms  in  the  system  are 
similar,  or  when  the  composition  includes  slightly  more  copper 
atoms.  This  suggests  that  good  GFA  in  this  system  is  a  balance 
between  having  as  many  atoms  of  each  of  the  three  elements  as 
possible,  and  having  a  large  number  of  atoms  of  copper,  which  has 
a  significantly  smaller  radius  than  calcium  and  magnesium  (see 
Table  1 ).  This  maximises  the  complexity  of  the  structure,  in  line 
with  the  empirical  rules  for  high  GFA  laid  down  by  Inoue  [17,18]. 

The  interpretation  of  diffraction  results  for  metallic  glass 
systems  which  contain  more  than  two  atom  types  is  challenging. 
The  atomic  radii  of  the  elements  involved  are  often  distributed 
such  that  the  first  peaks  of  the  partial  correlation  functions 
overlap  and  the  determination  of  individual  interatomic  distances 
from  the  total  correlation  function,  T(r),  is  extremely  difficult.  A 
common  approach  is  to  use  modelling  techniques,  such  as  RMC, 
where  models  are  constructed  to  fit  reciprocal  space  data,  and 
then  interatomic  distances  and  partial  coordination  numbers  are 
calculated  from  the  model.  A  recent  RMC  study  of  the  related 
Ca— Mg— Zn  glass  system  has  given  insights  into  the  changes  in 
glass  structure  with  composition,  and  how  this  affects  the  GFA, 
encouraging  crystallisation  [16].  The  direct  approach  for  extracting 
information  on  the  short  range  order  in  the  glass  is  to  fit  peaks  to 
the  correlation  function,  T(r),  but  for  ternary  glasses  there  are  six 
independent  partial  pair  correlation  functions  which,  in  principle, 
require  six  independent  experiments  with  different  elemental 
weighting  factors  for  a  unique  separation  of  partial  functions. 
Nonetheless,  the  fitting  method  has  been  applied  successfully  to 
the  ternary  glass  system  AlgoFexCeio-x  [19,20],  where  the 
compositions  are  favourable  for  the  fitting  technique  due  to  the 
low  concentration  of  iron  and  cerium  in  these  samples.  T(r)  is 
dominated  by  aluminium  correlations  and,  if  Fe  and  Ce  are 
considered  together  as  the  solute,  then  the  ECP  model  indicates 


Table  1 

The  metallic  (46]  and  covalent  (43]  radii  for  calcium,  magnesium  and  copper. 


Atom 

Rm(A) 

Rc(A) 

Ca 

1.97 

1.76 

Mg 

1.60 

1.41 

Cu 

1.28 

1.32 

that  there  will  be  no  Fe— Fe,  Ce— Ce  and  Fe— Ce  nearest  neighbours. 
Therefore,  the  fit  was  reduced  to  three  peaks  and  the  environ 
ments  of  all  three  elements  were  measured  reliably.  The  work 
carried  out  on  the  systems  Zr7oAlio(Ni,  Cu)2o  [21,22]  is  more 
relevant  because  the  concentration  of  solute  atoms  is  comparable 
to  the  glasses  studied  here.  The  multiple  peak  fitting  technique 
was  applied  to  the  X  ray  and  neutron  correlation  functions  with 
some  success  and  it  was  demonstrated,  with  additional  informa 
tion  from  Zr,  Ni,  and  Cu  EXAFS,  that  the  addition  of  aluminium  has 
a  large  impact  on  the  glass  structure  when  substituted  into  the 
binary  glasses  Zr7o(Ni,Cu)3o.  The  Zr  environment  was  determined 
for  both  glasses,  though  again  the  solute— solute  correlations  were 
neglected. 

In  this  study  the  structure  of  the  glass  system,  CagoMg^-xCux,  is 
examined  by  using  direct  fitting  of  total  correlation  functions,  T(r), 
measured  using  both  neutron  and  X  ray  diffraction.  For  this 
method  there  are  difficulties  which  arise  from  the  large  number  of 
overlapping  contributions,  and  the  fact  that  there  is  insufficient 
information  for  an  unambiguous  separation  of  partials.  However,  it 
is  demonstrated  that  important  structural  insights  can  be  gained  by 
using  this  method,  even  in  cases  where  solute— solute  correlations 
cannot  be  neglected,  provided  that  the  analysis  is  carried  out  in 
a  critical  manner,  and  with  reasonable  constraints  applied  to  the 
peak  positions  for  the  various  correlations.  To  assess  the  reliability 
of  the  coordination  numbers  extracted,  the  data  are  compared  to 
partial  correlation  functions,  tU'(r),  calculated  using  the  Per 
cus— Yevick  (PY)  approximation  for  a  binary  hard  sphere  system 
[23,24], 

2.  X-ray  and  neutron  diffraction  theory 

A  total  scattering  experiment  measures  the  total  (X  ray  or 
neutron)  scattering  from  the  sample,  /(Q),  where  Q  is  the  magni 
tude  of  the  scattering  vector  for  elastic  scattering  [25].  The  total 
scattering  is  the  sum  of  the  self  scattering,  f(Q)  (the  interference 
between  waves  scattered  from  the  same  atom),  and  the  distinct 
scattering,  i(Q)  (the  interference  between  waves  scattered  from 
different  atoms),  which  is  expressed  simply  for  neutron  diffraction 
as  a  differential  cross  section: 

Jg  iN(Q)  /s(Q)  +  *n(Q)  (i) 

For  both  X  rays  and  neutrons,  the  self  scattering  can  be  calcu 
lated  within  an  approximation  for  a  sample  of  known  composition, 
and  is  subtracted  from  the  data  to  give  the  distinct  scattering. 

For  neutron  diffraction,  iN(Q),  can  then  be  Fourier  transformed 
to  give  a  real  space  correlation  function, 

(\  2  00 

EC'M  +§/  QiN(Q)M(Q)sin(rQ)dQ  (2) 

'  '  o 

where  M(Q)  is  a  modification  function  which  is  used  to  reduce 
termination  ripples  in  the  Fourier  transform.  p°  is  the  atomic 
number  density,  and  c/  and  £>(  are  respectively  the  atomic  fraction 
and  coherent  neutron  scattering  length  for  element  /.  The  resolu 
tion  in  real  space  depends  on  the  modification  function,  and 
principally  on  the  maximum  momentum  transfer,  Qmax.  of  the 
experimental  data. 

In  the  case  of  X  ray  diffraction,  the  situation  is  more  complex 
because  X  rays  scatter  from  the  electron  cloud,  yielding  a  relatively 
broad  ‘electron— electron’  correlation  function.  To  account  for  this, 
the  sharpened  distinct  scattering  is  obtained  by  first  dividing  by 
(f(Q))2  (where  (f(Q))  is  the  mean  X  ray  form  factor  for  the  sample) 
[26,27],  to  approximate  scattering  from  point  sources; 
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ix(Q) 


JX(Q)  (f(Q))2 
(f(Q)}2 


(3) 


The  X  ray  correlation  function  is  then  obtained  according  to, 


Tx{r)  4np°r+2j  Qrx(Q)iW(Q)sin(rQ)dQ  (4) 

o 


The  contribution  to  the  distinct  scattering  due  to  a  single 
interatomic  distance,  r^,  with  a  root  mean  square  (RMS)  variation 
in  distance  of  (u^)1/2  is, 


njkcjbjbk 


sin  (Qrjk) 

Qrjk 


exp 


i « 

>q2\ 

l  2 

) 

n  jjmm  sin(Qr^) 
JkJ  m .))2  &jk 


exp 


(5) 


(6) 


were  collected  using  the  GEM  diffractometer  [31  ]  at  the  ISIS  pulsed 
neutron  source  at  the  Rutherford  Appleton  Laboratory,  UK,  and  the 
data  reduction,  correction  and  analysis  were  carried  out  using 
a  combination  of  the  program  GUDRUN  [32]  and  the  ATLAS  suite  of 
programs  [33]. 

The  X  ray  diffraction  data  were  collected  on  a  Panalytical  X’pert 
Pro  Multi  Purpose  Diffractometer  which  had  been  optimised  for 
the  study  of  disordered  materials  with  a  capillary  stage  for 
mounting  the  samples,  and  a  tube  with  a  silver  anode  for  short 
wavelength  X  rays  (0.560885  A).  The  samples  were  finely 
powdered  and  loaded  in  silica  capillaries  1  mm  in  diameter  and 
with  walls  of  a  nominal  thickness  of  0.01  mm.  The  data  were 
collected  using  a  silicon  scintillation  point  detector  on  a  spinning 
sample  for  scattering  angles  over  a  range  from  3.2°  to  128°  with 
0.2°  intervals.  The  data  were  corrected  for  background,  absorption, 
polarisation,  multiple  scattering,  Compton  scattering  and  brems 
strahlung  effects,  and  normalised  using  a  newly  developed 
program,  GudrunX  [34]. 

Simultaneous  fitting  of  the  X  ray  and  neutron  correlation  func 
tions,  Tx(r)  and  T^r)  was  carried  out  using  the  program  NXfit  [35], 


where  rijk  is  the  coordination  number  of  k  type  atoms  around  the 
atom  j,  and  Cj  is  the  atomic  fraction  for  j  type  atoms.  The  contri 
bution  to  the  correlation  function  due  to  the  interatomic  distance 
rjk  is  then  evaluated  by  Fourier  transformation  of  i^(Q)  or  ixk(Q), 
according  to  Equation  (2)  or  (4). 

The  coordination  numbers  for  j— k  and  k—j  atom  pairs  are  related 
by  the  following  identity  [28]: 

njkcj  n  kjck  (7) 

A  coordination  number  only  has  meaning  if  the  range  in  inter 
atomic  distances,  rjk,  over  which  it  applies,  and  the  method  by 
which  it  has  been  determined,  is  defined  [29,30).  Provided  that  rijk 
and  nkj  have  been  determined  over  the  same  distance  range,  and  by 
the  same  method,  Equation  (7)  is  always  true.  It  should  also  be 
noted  that  the  coordination  number  derived  from  the  correlation 
function  is  an  average  over  all  atomic  sites  in  the  sample. 

3.  Experimental  detail 

Three  samples  of  compositions,  Ca6oMgi5Cu25,  Ca6oMg2oCu2o 
and  Ca6oMg25Cui5  were  prepared  by  melt  spinning  and  charac 
terised  (including  density  measurements)  as  described  by  Senkov 
et  al.  [13], 

For  neutron  diffraction,  small  fragments  of  the  sample  were 
loaded  into  10.3  mm  diameter  cylindrical  vanadium  containers, 
with  walls  which  were  25  pm  in  thickness  to  reduce  background 
scattering  and  absorption  of  neutrons  by  the  container.  The  data 


4.  Results 

The  corrected  neutron  distinct  scattering,  iN(Q),  for  each  of  the 
three  samples  is  shown  in  Fig.  1(a),  with  expanded  plots  of  the  low 
Q  region  in  Fig.  2(a)  and  the  neutron  correlation  function,  T^r),  in 
Fig.  3(a).  For  each  figure,  the  X  ray  diffraction  data  are  shown  in 
part  (b),  adjacent  to  the  neutron  data  for  comparison.  ix(Q )  for  each 
of  the  three  samples  is  shown  in  both  Figs.  1(b)  and  2(b),  whilst 
Tx(r)  is  shown  in  Fig.  3(b).  The  correlation  functions  were  obtained 
by  Fourier  transformation,  according  to  Equations  (2)  and  (4),  using 
the  Lorch  modification  function  [36]  with  a  maximum  momentum 
transfer  Qmax  of  25  and  18  A-1  for  neutron  and  X  ray  diffraction 
respectively.  Differences  between  T^r)  and  Tx(r)  provide  more 
information  on  the  various  interatomic  correlations  in  the  glass 
than  can  be  obtained  by  the  use  of  one  radiation  alone. 

The  distinct  scattering,  i(Q),  for  both  radiations  damps  very 
quickly  with  Q  and  no  oscillations  are  apparent  beyond  15  A-1  (see 
Fig.  1 ),  indicating  that  the  glasses  have  a  wide  range  of  interatomic 
distances.  However,  small  Bragg  peaks  in  iN(Q )  (at  about  2.6  and 
3.7  A-1)  indicate  that  the  samples  are  not  completely  amorphous. 
There  is  a  small  amount  of  crystallisation  for  all  three  samples  and 
the  Bragg  peak  intensity  increases  as  copper  is  added.  The  peaks 
have  been  identified  as  arising  from  CaO  [37],  which  forms  on  the 
surface  of  the  glasses  (see  Fig.  2).  The  Bragg  peaks  are  observed 
more  clearly  in  the  neutron  diffraction  patterns  than  in  the  X  ray 
diffraction  patterns  because  the  neutron  diffraction  measurement 
has  a  narrower  resolution.  Fig.  2  shows  that  there  are  systematic 
changes  with  composition  in  the  distinct  scattering  which  are 


Fig.  1.  Distinct  scattering  data  for  three  glass  compositions  measured  using  a)  neutrons  and  b)  X-rays. 
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Fig.  2.  A  comparison  of  the  low-Q  region  of  i(Q)  for  the  three  samples.  Inset  is 
a  comparison  of  the  neutron  diffraction  data  for  Ca6oMgi5Cu25  and  the  Bragg  peaks 
expected  for  crystalline  CaO  [37], 

consistent  between  the  two  radiations.  As  copper  is  substituted 
into  the  glass  at  the  expense  of  magnesium,  the  first  sharp 
diffraction  peak  (FSDP)  shifts  to  higher  Q  and  broadens.  The 
correlation  functions  show  that  an  increase  in  the  Cu  content 
causes  a  shortening  of  the  average  interatomic  distance,  and  an 
increase  in  the  width  of  the  distribution  of  interatomic  distances. 
Thus  the  shift  of  the  FSDP  may  be  interpreted  as  arising  from  the 
shortening  of  the  average  interatomic  distance,  whilst  the  broad 
ening  of  the  FSDP  may  arise  from  the  broadening  of  the  distribution 
of  interatomic  distances. 

The  correlation  functions  show  a  peak  manifold  in  the  region 
2.2— 4.5  A,  which  arises  from  the  interatomic  distances  between  the 
various  atom  pairs  in  these  glasses.  Systematic  changes  with 
composition  are  also  observed  in  the  total  correlation  functions 
(see  Fig.  3).  As  the  amount  of  copper  in  the  sample  is  increased, 


there  is  an  increase  in  area  between  2.0  and  3.2  A,  and  a  corre 
sponding  decrease  in  area  between  3.2  and  4.0  A.  Fig.  4  shows  the 
differences  between  the  correlation  functions  with  the  two  most 
disparate  compositions  (Ca6oMgi5Cu25— Ca6oMg25Cuis).  The  posi 
tions,  of  the  maxima  and  minima  in  the  differences  have  been 
used  to  estimate  the  interatomic  distances  between  different 
elements  (see  Table  2).  For  both  radiations,  the  data  for  the  sample 
with  lower  Cu  content  are  subtracted  from  the  data  for  the  sample 
with  higher  Cu  content,  resulting  in  positive  peaks  arising  from 
copper  correlations,  and  negative  peaks  due  to  interatomic 
distances  involving  magnesium.  The  positions  of  the  first  three 
maxima  and  minima  can  be  determined  easily,  because  there  is 
good  agreement  between  the  results  from  the  two  radiations  and 
the  features  are  well  defined.  At  distances  greater  than  3  A, 
although  it  is  still  possible  to  make  an  assignment  of  the  inter 
atomic  distances,  it  is  more  difficult  because  the  peaks  are  less 
distinct,  with  more  variation  between  the  differences  measured  by 
the  two  radiations.  In  Table  2  the  bonding  for  each  pair  correlation 
is  categorised  as  metallic  or  covalent,  depending  whether  the  value 
of  rdiff  is  closer  to  the  sum  of  metallic  or  covalent  radii,  rM  or  r£. 

For  each  of  the  three  samples,  peak  fits  were  simultaneously 
made  to  both  T^r)  and  7*(r)  using  the  program  NXfit  [35].  The 
presence  of  three  atom  types  with  similar  atomic  radii  (see  Table  1 ) 
results  in  six  different  correlations  with  nearest  neighbour 
distances  between  2  and  4  A  (see  Table  2),  and  it  is  impossible  to 
deconvolute  these  six  contributions  by  a  unique  fit  to  T(r).  However, 
by  utilising  the  different  scattering  powers  of  the  elements  for  the 
two  radiations,  and  fitting  the  two  correlation  functions  simulta 
neously,  the  range  of  possible  fits  to  the  data  is  drastically  reduced. 
When  this  method  of  fitting  is  combined  with  the  information  on 
interatomic  distances  derived  from  the  results  in  Fig.  4,  a  final  fit 
can  be  achieved  which  is  a  reasonable  model  for  the  short  range 
order  of  the  glass. 

The  initial  simultaneous  fit  was  made  to  the  correlation  func 
tions  for  Ca6oMg15Cu25.  Five  peaks  were  used  to  model  each  of  the 
possible  atom  pairs,  excluding  Mg— Mg.  Magnesium  weakly  scatters 
both  X  rays  and  neutrons  and  this,  coupled  with  the  small  number 
of  magnesium  atoms  present  in  the  three  samples,  leads  to 
a  contribution  to  T(r)  from  Mg— Mg  distances  which  can  be 
neglected  without  much  effect  on  the  measured  coordination 
numbers  for  the  other  correlations.  The  initial  positions  for  the  five 
peaks  were  derived  from  the  differences  shown  in  Fig.  4,  and  were 
then  allowed  to  vary  within  limited  bounds  (±0.05  A),  whilst  the 
peak  widths  and  areas  were  allowed  to  vary  without  constraint.  The 
resultant  fit  was  then  used  to  provide  the  initial  fit  parameters  for 
the  fit  to  the  Ca6oMg2oCu2o  data,  and  only  the  peak  area  (i.e.  coor 
dination  number)  was  allowed  to  vary.  This  resulted  in  a  poor  fit  in 


Fig.  3.  The  total  correlation  function  measured  with  neutrons  (a)  and  X-rays  (b).  The  differences  are  also  shown,  where  the  dashed  grey  line  is  the  difference  Ca60Mg15 
Cu25  Ca6oMg2oCu2o  and  the  grey  dotted  line  is  the  difference  Ca6oMgi5Cu25  Ca6oMg25Cui5. 
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Fig.  4.  Differences  between  the  total  correlation  functions  of  the  two  most  disparate 
compositions  (Ca6oMgi5Cu25  Ca6oMg25Cui5)  for  neutrons  (grey)  and  X-rays  (black). 
The  neutron  difference  has  been  multiplied  by  2.5  to  be  comparable  with  the  X-ray 
data. 

the  region  of  the  Mg— Ca  and  Ca— Ca  correlations,  and  so  the  peak 
positions  and  widths  of  those  peaks  were  then  also  allowed  to  vary. 
The  final  fitting  parameters  for  this  sample  were  then  used  in  the 
same  way  as  the  starting  parameters  for  the  fit  to  the  Ca6oMg25Cui5 
data.  Figs.  5—7  show  the  final  fits  to  the  X  ray  and  neutron  data  for 
each  of  the  samples. 

Once  satisfactory  fits  to  the  X  ray  and  neutron  datasets  were 
achieved,  the  coordination  numbers  (and  peak  widths,  where 
appropriate)  for  each  sample  were  altered,  and  the  program  was 
allowed  to  refit  the  data.  The  differences  between  the  final  values 
for  the  two  fits  were  then  used  as  an  estimate  of  the  error.  The  fit 
parameters  for  the  three  correlations  involving  copper  are  given  in 
Table  3.  The  Mg— Cu  and  Ca— Cu  coordination  numbers  were 
calculated  according  to  Equation  (7).  Values  for  the  Mg— Ca  and 
Ca— Ca  peaks  are  not  reported  as  they  are  unreliable  and  only 
present  to  improve  the  reliability  of  the  fits  to  the  copper  peaks. 
Fig.  8  shows  the  variation  with  composition  of  the  fitted  coordi 
nation  numbers  for  each  of  the  copper  correlations.  Each  coordi 
nation  number  was  found  to  vary  linearly  with  copper  content  over 
the  studied  composition  range. 

5.  Discussion 

5.3.  Peak  fits  and  coordination  numbers 

There  are  suggestions,  from  both  modelling  studies  [5,9]  and 
experimental  studies  [21,22],  that  there  may  be  no  solute— solute 
contacts  for  ternary  BMGs  with  low  solute  concentrations.  The  ECP 
model  allows  the  number  of  structural  sites  to  be  counted  [9],  and 

Table  2 

The  closest  approach  of  each  interatomic  pair,  rj?  and  rfk,  calculated  using  metallic 
[46]  and  covalent  [43]  radii  respectively.  These  distances  are  compared  to  the 
interatomic  distances,  rdiff,  determined  from  the  composition  differences  shown  in 
Fig.  4,  and  the  type  of  bonding  for  each  correlation  has  been  identified. 


Atom  pair,  j-k 

r“  (A) 

rfk  (A) 

rdiff  (A)  (Fig.  4) 

Bonding 

Cu  Cu 

2.56 

2.64 

2.5 

Metallic 

Cu  Mg 

2.88 

2.73 

2.7 

Covalent 

Cu  Ca 

3.25 

3.08 

3.0 

Covalent 

Mg  Mg 

3.2 

2.82 

Mg  Ca 

3.57 

3.17 

3.6 

Metallic 

Ca  Ca 

3.94 

3.52 

3.8 

Metallic 

r 


Fig.  5.  Fits  to  X-ray  and  neutron  correlation  functions  for  Ca60Mg15Cu25  (the  black  line 
shows  the  experimental  data,  whilst  the  dashed  line  shows  the  fitted  function,  and  the 
grey  lines  show  the  components  of  the  fit). 


so  gives  a  formalism  for  calculating  the  maximum  solute  concen 
tration  where  solute— solute  contact  can  be  avoided.  Consider 
a  Ca— Mg— Cu  glass  structure  where  all  solute  sites  are  occupied 
only  by  solute  atoms  and  each  solute  site  is  surrounded  only  by 
solvent  atoms.  Using  Mg  as  the  primary  structure  forming  solute 
(a),  the  relative  size  of  Mg  and  Ca  atoms  suggests  that  each  a  site 
will  be  surrounded  by  ~  10  Ca  atoms.  Three  additional  solute  sites 
occur  at  the  interstices  of  these  structure  forming  a  centred  clus 
ters.  There  are  thus  a  total  of  14  structural  sites  per  a  site,  and  4  of 
these  are  solute  sites.  The  maximum  solute  atom  fraction  where 
solute— solute  contact  can  be  avoided  is  thus  4/14  =  0.29.  If  a  sites 
are  primarily  occupied  by  Cu  atoms,  then  the  number  of  Ca  atoms 
surrounding  each  a  site  is  —  9  and  the  critical  solute  concentration 
below  which  solute— solute  contact  may  be  avoided  is  4/13  =  0.31. 
Solute— solute  contact  becomes  topologically  necessary  for  solute 
concentrations  greater  than  these  values. 

In  agreement  with  this  simple  analysis,  the  correlation  functions 
measured  for  the  calcium— magnesium— copper  system,  shown  in 
Fig.  3,  indicate  that  the  glass  structure  in  this  system,  which  is 
comprised  of  40%  solute  atoms,  can  only  be  explained  if  sol 
ute— solute  contact  occurs  in  these  glasses.  The  complexity  of  the 
peak  manifold  between  2.2  and  4.5  A,  and  the  changes  which  occur 
on  altering  the  relative  amount  of  copper  and  magnesium  with 
respect  to  calcium,  indicate  that  there  are  more  contributions  than 
can  arise  from  Cu— Ca,  Mg— Ca,  and  Ca— Ca  correlations  alone.  Thus, 
the  solute— solute  correlations  Cu— Cu,  Cu— Mg  and  Mg— Mg  cannot 
be  neglected.  While  glasses  in  the  Ca— Mg— Cu  system  can  be 
formed  with  compositions  similar  to  Zr7oAlio(Ni,  Cu)2o  (i.e. 
Ca7o(Mg,Cu)3o),  the  GFA  is  poor  [13].  This  suggests  that  higher 
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Fig.  6.  Fits  to  X-ray  and  neutron  correlation  functions  for  Ca6oMg2oCu2o  (the  black  line 
shows  the  experimental  data,  whilst  the  dashed  line  shows  the  fitted  function,  and  the 
grey  lines  show  the  components  of  the  fit). 


levels  of  solute,  and  solute— solute  correlations,  are  beneficial  to  the 
formation  of  a  BMG.  This  beneficial  effect  has  been  proposed  to 
come  from  an  increase  in  the  more  stable  solute— solvent  bonds 
that  accompany  an  increase  in  the  solute  content  [9], 

The  simultaneous  fitting  of  the  X  ray  and  neutron  correlation 
functions  achieves  good  agreement  with  the  experimental  data  by 
using  the  interatomic  distances  identified  from  Fig.  4  as  a  starting 
point.  The  narrow  widths  of  the  resolution  broadened  Gaussian 
peaks  modelling  the  Cu— Cu,  Cu— Mg  and  Cu— Ca  correlations  (see 
Figs.  5—7),  are  consistent  with  strong  bonding  of  the  copper.  In 
contrast,  the  two  fitted  peaks  for  the  Mg— Ca  and  Ca— Ca  correla 
tions  are  very  broad  and  it  was  necessary  to  allow  both  the  position 
and  width  of  the  peaks  to  vary  between  samples  to  achieve  a  better 
fit.  Although  broader  peaks  can  be  an  indication  of  increased  static 
disorder,  it  is  likely  that  in  this  case  the  fitted  peaks  are  not  accurate 
representations  of  the  distribution  of  Mg— Ca  and  Ca— Ca  inter 
atomic  distances.  Furthermore,  extra  area  in  the  region  of  3.2—4  A 
may  arise  from  a  second  shell  of  copper  correlations  (as  in  crys 
talline  Ca2Cu  —  see  below),  and  therefore  any  attempts  to  extract 
structural  information  from  the  fits  to  the  Mg— Ca  and  Ca— Ca  peaks 
are  less  likely  to  yield  reliable  information.  The  consistency  of  the 
fits  to  the  copper  correlations  indicates  that  realistic  coordination 
numbers  for  copper  can  be  measured  over  the  limited  range  of 
~2— 3.2  A  and  the  fitting  parameters  for  the  copper  related 
correlations  are  given  in  Table  3. 

For  a  glass  with  a  majority  of  calcium  atoms,  the  coordination 
number  for  a  minority  atom  can  be  calculated  using  both  the  dense 
random  packing  (DRP)  model  [38,39]  and  the  ECP  model  [5,9].  For 
a  calcium  based  glass  with  minority  copper  atoms,  both  models 
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Fig.  7.  Fits  to  X-ray  and  neutron  correlation  functions  for  Ca6oMg25Cui5  (the  black  line 
shows  the  experimental  data,  whilst  the  dashed  line  shows  the  fitted  function,  and  the 
grey  lines  show  the  components  of  the  fit). 

indicate  that  the  number  of  calcium  nearest  neighbours  to  a  copper, 
ticuca.  is  between  8  and  9.  However,  the  total  copper  coordination 
number  determined  from  the  fits  varies  from  4.74  for  the  sample 
with  the  lowest  copper  content,  to  5.74  for  the  sample  containing 
the  most  copper;  clearly  this  is  significantly  less  than  predicted  by 
the  DRP  or  ECP  models. 

The  crystal  structure  of  Ca2Cu  [40]  can  be  considered  as  a  binary 
crystalline  analogue  for  the  glass  because  it  has  a  similar  solute/ 
solvent  atom  fraction.  In  this  crystal  a  copper  atom  has  9  nearest 
neighbours  closer  than  3.5  A  (the  longest  Cu— Ca  distance  identified 
by  the  fits  shown  in  Figs.  5—7),  agreeing  closely  with  the  predictions 
of  the  DRP  [38,39]  and  ECP  models  [5,9],  though  these  models 
assume  that  all  the  neighbours  are  calcium  while  in  the  crystal 


Table  3 

The  fitted  parameters  for  the  copper  correlations  in  the  three  glasses.  The  inter¬ 
atomic  distances  were  fixed  for  all  three  samples.  The  peak  widths  were  allowed  to 
vary  for  the  Ca6oMgi5Cu25  sample,  and  were  subsequently  fixed.  Values  given  in 
square  brackets  and  Italics  are  calculated  from  the  PY  simulation  for  Ca60Cu4o. 


Atomic 
pair,  j-k 

Interatomic 

distance, 

(A) 

Peak 

width, 

K)V2 

(A) 

Ca6oMg15Cu25 
Coordination 
number,  njk 

Ca6oMg2oCu2o 
Coordination 
number,  njfe 

CaeoMg25Cui5 

Coordination 
number,  njk 

Cu  Cu 

2.52* 

0.12* 

1.18(5) 

0.88(6) 

0.74(6) 

Cu  Mg 

2.73* 

0.18* 

0.51(8) 

0.77(4) 

0.92(5) 

Cu  Ca 

3.04* 

0.15* 

4.05(4)  (4.85] 

3.56(6)  [4.85] 

3.08(6)  [4.85] 

Mg  Cu 

2.73* 

0.18* 

0.9(1) 

0.77(4) 

0.55(3) 

Ca  Cu 

3.04* 

0.15* 

1.69(7)  [2.02] 

1.19(6)  [1.62] 

0.77(3)  [1.21] 

Total  Cu  coordination  number 

5.74(9)  [6.54] 

5.21(8)  [6.50] 

4.74(8)  f 6.5 1 7 

Density  (g  cm  3) 

2.367(1) 

2.200(1) 

2.039(1) 
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Fig.  8.  Cu  coordination  numbers,  with  closed  symbols  indicating  Cu  X,  and  open 
symbols  indicating  X  Cu.  Squares  are  Cu  Cu,  circles  are  Cu  Mg  (and  Mg  Cu)  and 
triangles  are  Cu  Ca  (and  Ca  Cu).  The  continuous  lines  are  fits  of  the  form 
n jk  constant  x  q(,  whilst  the  dashed  lines  are  linear  fits.  The  symbols  (labelled  PY)  at 
40  at.%  Cu  indicate  the  coordination  numbers  obtained  by  fitting  resolution-broadened 
Gaussians  to  the  Cu  Ca,  Ca  Cu  and  Cu  Cu  partial  correlation  functions  calculated  for 
the  binary  glass  Ca6oCu40  using  the  Percus  Yevick  hard  sphere  approximation. 


there  are  two  Cu  neighbours,  as  shown  in  Fig.  9.  Furthermore,  the 
average  Cu— Ca  distance  is  3.07  A,  which  is  consistent  with  covalent 
bond  lengths,  and  is  much  shorter  than  3.24  A,  the  average  of  the 
Cu— Cu  (2.51  A)  and  Ca— Ca  (3.98  A)  distances.  This  is  in  good 
agreement  with  the  interatomic  distances,  r^ff  (Table  2),  observed 
in  Fig.  4  and  the  non  additive  nature  of  bond  lengths  in  BMGs  has 
been  discussed  previously  in  a  number  of  papers  including  one  by 
Cheng  et  al.  [41  ].  Fig.  9  shows  a  comparison  of  the  neutron  corre 
lation  function  for  the  glass  (Ca6oMgi5Cu25)  of  closest  composition 
to  the  crystal,  with  a  simulation  of  the  three  contributions,  Cu— Cu, 


r 

Fig.  9.  A  comparison  the  experimental  7 w(r)  for  Ca6oMg15Cu25  with  the  Cu  Cu,  Cu  Ca 
and  Ca  Ca  contributions  simulated  (see  text  for  details)  according  to  the  Ca2Cu  crystal 
structure  [40],  The  inset  shows  the  atomic  arrangement  around  a  Cu  atom  in  the  crystal, 
including  all  atoms  closer  than  3.5  A  (small  spheres  are  Cu,  large  spheres  are  Ca). 


Cu— Ca  and  Ca— Ca,  based  on  the  Ca2Cu  crystal  structure  [40],  For  the 
simulation,  the  partial  correlation  functions  for  the  crystal  were 
calculated  using  the  XTAL  program  [42],  and  were  then  weighted 
appropriately  for  the  composition  of  the  glass.  The  peak  width  used 
to  broaden  each  of  the  partials  was  0.15  A,  a  typical  value  for  the  fits 
to  the  experimental  data  given  in  Table  3.  The  positions  of  the 
simulated  Cu— Cu  and  Cu— Ca  peaks  are  similar  to  those  for  the  glass, 
but  the  peaks  are  too  intense,  indicating  that  the  coordination 
numbers  for  Cu— Cu  and  Cu— Ca  in  the  glass  are  significantly  less 
than  2  and  7  respectively,  in  agreement  with  the  fit  results  in  Table 
3. 

It  should  be  noted  that  the  covalent  Cu— Cu  distance  in  Table  2  is 
unusual  in  that  it  is  longer  than  the  metallic  distance.  The  tabulated 
covalent  Cu  radius  was  calculated  as  an  average  from  a  survey  of 
many  crystal  structures  [43],  This  unusual  behaviour  may  be 
related  to  the  electronic  structure  of  a  copper  atom,  which  has 
a  single  electron  in  the  outermost  orbital.  This  electron  is  donated 
to  form  metallic  bonds  and  reduces  the  atomic  radius  to  a  greater 
extent  than  occurs  when  atoms  overlap  to  form  covalent  bonds.  The 
good  agreement  between  the  Cu— Mg  and  Cu— Ca  distances  calcu 
lated  using  the  literature  values  for  covalent  radii  and  the  distances 
determined  from  the  diffraction  data  gives  confidence  that  the 
values  used  for  the  covalent  radii  (notably  for  Cu)  are  reasonable. 

The  three  Cu— X  coordination  numbers  for  each  sample  (and  the 
corresponding  X— Cu  coordination  numbers,  calculated  according 
to  Equation  (7))  are  shown  as  a  function  of  copper  content  in  Fig.  8. 
For  all  five  coordination  numbers,  the  changes  with  composition 
are  significantly  larger  than  the  errors  from  varying  the  fitting  of 
the  data  (indicated  by  error  bars  in  Fig.  8).  Within  errors,  the  sol 
ute— solute  coordination  numbers,  are  proportional  to  the  relevant 
atomic  fraction,  so  that  ricucu  (4.7  ±  0.3)ccu,  riMgCu  (3.6  ±  0.2) 
ecu,  and  nCUMg  (3.6  ±  0.2)cMg,  as  indicated  by  the  continuous  lines 
in  Fig.  8.  Cargill  &  Spaepen  [28]  have  stated  correctly  that  the 
coordination  number  of  a  system  with  no  chemical  order  is  only 
proportional  to  if  the  two  atom  types,  j  and  k,  are  the  same  size, 
and  hence  our  results  initially  seem  surprising.  To  further  test  these 
results,  Monte  Carlo  simulations  for  model  Ca6oMg4o_xCux  glasses 
with  a  packing  fraction  of  50%  and  hard  sphere  potentials  based  on 
the  metallic  radii  given  in  Table  1  were  performed  [44],  These 
models  show  that  for  no  chemical  order,  the  coordination  numbers, 
rijh,  where  k  is  a  solute  element,  are  proportional  to  c/<.  This  may  be 
interpreted  as  arising  because  the  solute  atoms  fill  holes  in  the 
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solvent  atom  matrix,  and  these  holes  are  the  same  size  irrespective 
of  the  solute  atom  type  and  able  to  accept  either  magnesium  or 
copper  atoms  without  causing  a  significant  restructuring  of  the 
atomic  packing.  In  this  case,  it  may  be  appropriate  to  consider  these 
three  glasses  as  a  pseudo  binary  glass  Ca6o(Cu/Mg)4o. 

For  the  Ca— Cu  coordination  numbers,  the  continuous  line  in 
Fig.  8  indicates  a  fit  based  on  the  assumption  of  proportionality  to 
the  Cu  atom  fraction,  ccu;  clearly  the  values  of  ncacu  determined 
from  fitting  the  experimental  data  are  not  proportional  to  ccu  (if  the 
values  of  ncacu  were  proportional  to  ccu,  then  Equation  (7)  shows 
that  McuCa  would  be  constant,  as  indicated  by  a  continuous  line  in 
Fig.  8).  The  dashed  lines  in  Fig.  8  are  linear  fits  to  the  data  for  the 
atom  pair  Ca— Cu,  with  the  corresponding  equations  in  Table  4. 
Fiowever,  these  linear  fits  are  not  satisfactory,  firstly  because  the  fit 
to  ncacu  becomes  negative  for  low  Cu  content,  and  secondly  because 
the  two  fits  are  fundamentally  incorrect  in  that  they  do  not  obey 
Equation  (7).  Furthermore,  the  values  obtained  for  nCuca  do  not 
behave  as  expected;  the  argument  that  the  three  glass  composi 
tions  can  be  considered  as  a  pseudo  binary,  necessitates  that  direct 
substitution  of  copper  for  magnesium  atoms  should  have  little  or 
no  effect  on  the  packing  of  calcium  atoms  around  copper.  Fiowever, 
the  fitting  results  indicate  that,  as  copper  atoms  are  added  to  the 
glass,  n cuca  increases.  The  fit  to  the  Ca— Cu  peak  is  the  least  reliable 
of  the  three  copper  peaks,  because  it  is  at  the  longest  distance,  and 
is  more  affected  by  overlap  problems  with  longer  distance  peaks, 
such  as  the  ‘dummy’  Mg— Ca  peak.  Therefore  the  values  given  by 
peak  fitting  should  be  treated  with  caution.  However,  further 
insight  into  the  Ca— Cu  peak  may  be  obtained  by  consideration  of 
hard  sphere  structural  models,  as  discussed  below. 

5.2.  Hard  sphere  calculations 

5.2.1.  Modelling  T(r)  using  the  Percus—Yevick  approximation 

In  order  to  gain  further  insight  into  the  coordination  numbers 
and  distributions  of  interatomic  distances,  and  their  derivation 
from  experimental  diffraction  data,  partial  correlation  functions 
were  calculated  for  the  hypothetical  binary  hard  sphere  systems, 
Ca6oCu4o  and  Ca6oMg4o  using  the  Percus—Yevick  (PY)  approxima 
tion  [23,24].  A  PY  simulation  of  a  hard  sphere  system  is  considered 
to  be  an  appropriate  model  for  the  structure  of  this  glass  system 
because,  by  definition,  hard  spheres  have  no  chemical  interactions 
and  hence  the  packing  is  random,  in  agreement  with  the  sol 
ute— solute  peak  fitting  results.  For  this  consideration,  the  Simula 
tions  are  compared  with  only  the  neutron  diffraction  data,  because 
of  the  superior  real  space  resolution  for  the  neutron  diffraction 
data. 

The  PY  calculations  directly  yield  partial  structure  factors,  and 
these  were  multiplied  by  a  Debye— Waller  factor,  exp(  u2/2Q2),  to 
represent  the  effect  of  atomic  vibrations,  and  then  Fourier  trans 
formed  (as  in  Equation  (2)),  using  exactly  the  same  modification 
function  as  for  the  experimental  data,  to  obtain  the  partial  corre 
lation  functions,  tfj!(r),  which  describe  correlations  between  an 
average  origin  atom  of  element  i,  and  atoms  of  element  /'  at 
a  distance  r.  For  the  PY  calculations  it  is  necessary  to  specify  the 
composition,  the  hard  sphere  radii  and  the  density  (or  packing 
fraction),  in  order  to  perform  calculations  for  a  hypothetical 
Ca6oCu40  glass,  the  density  was  estimated  as  2.858  g  cnrT3  by  linear 


Table  4 

Linear  fits  to  the  dependence  on  Cu  content  of  coordination  number  data  (Fig.  8). 


Atomic  pair 

Gradient 

Intercept 

Ca  Cu 

0.090  ±  0.007 

-0.59  ±  0.12 

Cu  Ca 

0.097  ±  0.007 

1.62  ±  0.15 

extrapolation  of  the  ternary  densities  (Table  3),  leading  to  a  packing 
fraction  of  0.61  which  was  used  for  all  the  PY  calculations.  The 
pairwise  sums  of  the  hard  sphere  radii  give  the  absolute  lower  limit 
of  approach  for  the  hard  spheres,  whereas  the  atomic  radii  given  in 
Table  1  are  related  to  the  mean  interatomic  distances.  Hence  it  is 
not  appropriate  to  use  the  atomic  radii  as  the  hard  sphere  radii  in 
the  PY  calculations.  It  was  found  empirically  that  hard  sphere  radii 
that  are  6.25%  smaller  than  those  given  in  Table  1  give  peak  posi 
tions  which  both  approximate  the  interatomic  distances  in  Table  2, 
and  simulate  the  peaks  in  the  experimental  correlation  functions 
well.  An  RMS  thermal  variation  in  interatomic  distance,  u,  of  0.1  A 
was  found  to  give  similar  peak  widths  to  the  experimental  data.  The 
values  in  Table  3  are  somewhat  larger  than  this  value,  since  they 
represent  the  combined  effects  of  thermal  and  static  disorder. 

Fig.  10  shows  a  comparison  of  the  experimental  neutron  corre 
lation  functions,  T*{r),  with  a  simulation  based  on  the  PY 


r 

Fig.  10.  A  comparison  of  T N(r)  for  the  three  glasses  with  a  simulation  based  on  the  hard 
sphere  PY  calculations  for  Ca6oCu4o  (see  text  for  details). 
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calculations  for  a  hard  sphere  Ca60Cu4o  glass.  The  simulated  T^r) 
was  calculated  from  a  sum  of  the  partial  correlation  functions  in 
which  they  were  weighted  according  to  the  composition  of  each 
ternary  glass.  This  makes  the  assumption  that  the  coordination 
numbers  for  Cu— Cu,  Cu— Ca  and  Ca— Ca  are  unchanged  as  Mg  is 
substituted  for  Cu.  In  the  regions  of  the  three  peaks  (between  2.3 
and  3.0  A  and  around  4.0  A),  the  simulation  gives  a  reasonably  close 
description  of  the  experimental  results,  but  the  lack  of  correlations 
involving  magnesium  in  the  model  leads  to  a  lack  of  intensity  in  the 
region  between  3  and  4  A  where  Mg— Mg  and  Mg— Ca  correlations 
are  expected. 

As  indicated  in  Table  2,  the  Cu— X  (whereX  is  not  Cu)  distances  in 
these  glasses  are  consistent  with  the  sum  of  covalent  radii,  whereas 
all  other  distances  are  consistent  with  the  sum  of  metallic  radii. 
Thus  Fig.  11(a)  shows  the  PY  simulation  of  the  Cu— Cu  partial 
correlation  function,  generated  using  metallic  radii,  whilst  Fig.  11(b) 
shows  the  PY  simulation  of  the  Ca— Cu  (or  Cu— Ca)  partial  correla 
tion  function  calculated  using  covalent  radii.  It  should  be 
acknowledged  that  the  simulation  of  different  partial  correlation 
functions  using  different  sets  of  radii  does  not  correspond  to  a  self 
consistent  model.  However,  this  approach  does  provide  a  reason 
able  estimate  of  the  form  of  the  distribution  of  nearest  neighbour 
interatomic  distances  for  each  pair  of  elements,  and  it  also  gives  an 
estimate  of  longer  range  correlations.  In  particular,  the  first  peak  of 
each  simulated  partial  correlation  function  is  asymmetric  with 
a  sharper  leading  edge,  and  a  broader  trailing  edge.  This  peak  shape 
is  consistent  with  some  real  space,  high  resolution  experimental 
results  on  the  partial  correlation  functions  for  metallic  glasses  (for 
example,  see  Ref.  [45]).  In  terms  of  a  hard  sphere  model,  the  sharp 
leading  edge  can  be  interpreted  as  arising  from  spheres  which  are 
in  direct  contact,  whilst  the  broad  trailing  edge  arises  from  spheres 
which  are  almost  touching,  but  are  not  actually  in  contact  because 
their  positions  are  constrained  by  all  the  other  spheres  which  are 
packed  around  them.  The  fitting  to  the  diffraction  data  of  a  single 
symmetric  peak  for  each  atom  pair  corresponds  more  closely  to  the 
contribution  to  the  correlation  function  from  atoms  which  are  in 
direct  contact.  However,  the  traditional  concept  of  a  ‘first  coordi 
nation  shell’  may  perhaps  correspond  more  closely  to  the  combined 
contributions  from  atoms  which  are  almost  touching,  as  well  as 
those  which  are  in  direct  contact  (Fig.  12). 

5.2.2.  Extracting  coordination  numbers  from  the  Percus—Yevick 
simulation 

The  coordination  numbers  for  the  PY  simulation  of  Ca6oCu4o 
were  evaluated  in  two  ways:  Firstly,  for  each  atom  pair,  a  resolu 
tion  broadened  Gaussian  was  fitted  to  the  leading  edge  of  the  first 


peak  in  tj)J(r)  (as  shown  in  Fig.  10),  and  secondly  rt[f,(r )  was  inte 
grated  up  to  the  minimum  after  the  first  peak.  The  first  approach 
gives  a  closer  simulation  of  the  way  in  which  the  experimental  data 
were  fitted,  and  only  includes  neighbours  which  are  in  close 
contact,  whereas  the  second  approach  gives  the  total  coordination 
number  for  the  first  coordination  shell,  for  distances  up  to  the 
minimum  in  the  distribution.  The  coordination  numbers  for  both 
approaches  are  given  in  Table  5. 

The  value  of  nCucu  from  the  fit  to  the  PY  simulation  agrees, 
within  error  (see  Table  5),  with  the  coordination  number  from 
linear  extrapolation  of  the  experimental  coordination  numbers, 
using  the  proportionality  relation  (continuous  line  in  Fig.  8).  This  is 
a  further  indication  that  the  Cu— Cu  peaks  fitted  to  the  glass  data  are 
an  accurate  representation  of  the  close  contact  part  of  the  distri 
bution  of  Cu— Cu  distances  for  randomly  distributed  copper  atoms. 
The  solute— solute  coordination  numbers,  n jk,  derived  by  fitting  the 
experimental  data,  have  a  physically  reasonable  behaviour,  in  that 
they  are  proportional  to  Ck,  and  the  agreement  between  the 
extrapolation  of  the  fit  values  for  ncucu  and  the  PY  result  is  further 
confirmation  that  the  values  are  physically  reasonable.  The  Cu— Cu 
component  of  the  fits  to  the  experimental  data  is  the  most  reliable 
of  all  components  because  it  occurs  at  shortest  distance,  and  this  is 
why  it  yields  reasonable,  reliable  results. 

As  noted  in  the  previous  section,  the  behaviour  of  the  Ca— Cu 
coordination  numbers  obtained  from  fitting  the  data  is  not  physi 
cally  reasonable.  One  of  the  reasons  for  this  is  that  the  fit  to  the 
Ca— Cu  peak  in  the  experimental  data  is  strongly  affected  by  the 
relatively  broad  ‘dummy’  Mg— Ca  peak  (see  Figs.  5—7).  Thus  it 
would  not  be  justifiable  to  extrapolate  the  values  for  the  Ca— Cu 
coordination  numbers  in  the  same  way  as  for  the  Cu— Cu  coordi 
nation  numbers,  and  instead  a  different  approach  was  adopted. 
This  approach  is  to  start  with  the  value  of  3.23  for  the  Ca— Cu 
coordination  number  from  the  PY  simulation  for  Ca6oCu40,  and  to 
extrapolate  to  the  compositions  of  the  glass  samples  on  the 
assumption  that  all  three  glasses  have  a  very  similar  arrangement 
of  the  solvent  atom  matrix  as  a  pseudo  binary  Ca6o(Ca/Mg)4o.  The 
results  from  this  approach  are  shown  in  Fig.  12  and  the  discrepancy 
between  them  and  the  experimental  coordination  numbers  grows 
as  the  Mg  content  increases,  due  to  the  increasing  effect  from 
overlap  with  the  Mg— Ca  peak  (see  Figs.  5—7).  Therefore,  by  using 
the  PY  simulation,  more  reasonable  values  for  nCacu  and  ncuca  have 
been  estimated  than  were  achievable  by  peak  fitting.  A  comparison 
of  the  coordination  numbers  is  given  in  Table  3. 

The  values  for  the  total  copper  coordination  numbers  calculated 
using  the  PY  derived  values  for  nCuca  and  nCa cu  are  given  in  Table  3 
(square  brackets),  and  are  significantly  larger  than  measured  by 


Fig.  11.  Partial  correlation  functions,  tjj[(r),  calculated  using  the  Percus  Yevick  approximation  for  a  binary  hard  sphere  system,  showing  a)  the  Cu  Cu  function,  and  b)  the  Cu  Ca 
function,  with  fits  to  the  leading  edge  of  a  resolution-broadened  Gaussian. 
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Fig.  12.  nCacu  and  nCucj  predicted  on  the  basis  of  the  PY  simulation  for  Ca60Cu40  and 
assuming  random  packing  with  no  chemical  effects  (circles),  together  with  the  results 
from  fitting  the  experimental  data  (triangles).  The  lines  show  the  composition- 
dependence  of  the  coordination  numbers  according  to  this  approach. 

peak  fitting  and  do  not  vary  significantly  with  composition.  The 
values  are  still  less  than  expected  from  the  DRP  and  ECP  models,  but 
these  results  are  a  measurement  of  the  symmetric  distribution  of 
the  nearest  neighbours  only.  The  PY  simulations  indicate  that  the 
distributions  of  interatomic  distances  are  likely  to  be  asymmetric, 
and  the  coordination  numbers  involving  copper  (Table  3)  are  likely 
to  be  an  underestimate  of  the  true  coordination  numbers,  due  to 
the  failure  to  include  the  broader  trailing  edge  in  the  peak  fitting. 
Table  5  shows  that  integrating  the  total  area  under  the  first  peak  in 
tcuCa(f).  calculated  by  PY,  gives  a  value  for  ncuca  of  8.83,  which  is 
much  closer  to  the  coordination  number  expected  from  consider 
ation  of  the  DRP  and  ECP  models  for  binary  glasses,  and  from  the 
crystal  Ca2Cu. 

From  the  work  presented  here  (and  elsewhere  [45])  it  appears 
that  a  PY  hard  sphere  simulation  may  provide  a  realistic  prediction 

Table  5 

The  Cu  Cu,  Cu  Ca  and  Ca  Cu  coordination  numbers  for  Ca60Cu4o,  obtained  by 
linear  extrapolation  of  the  results  from  peak  fitting  the  experimental  data  for  ternary 
glasses.  These  are  compared  with  values  obtained  from  PY  binary  hard  sphere 
simulations  of  the  partial  correlation  functions  of  Ca60Cu40  by  peak  fitting  and  by 
integration  (see  text  for  details). 


Atomic 
pair,  j  k 

njk  from  extrapolation  of 
parameters  from  peak  fits 
to  experimental  data  (with  errors) 

rijk  from 
PY 

peak  fit 

njk  from 

PY 

integration 

Cu  Cu 

1.87  (12) 

1.77 

3.82 

Cu  Ca 

5.50  (32) 

4.85 

8.83 

Ca  Cu 

3.05  (30) 

3.23 

5.89 

Cu  Cu  +  Cu  Ca 

7.37 

6.62 

12.65 

of  the  coordination  numbers,  and  notably  the  distribution  of 
interatomic  distances,  in  metallic  glasses  without  significant 
chemical  ordering.  However,  our  results  concerning  metallic  and 
covalent  radii  (Table  2)  show  that  it  would  be  advantageous  to  have 
a  theoretical  solution  of  the  PY  equation  for  a  binary  or  ternary  hard 
sphere  system  in  which  the  interatomic  distances  for  the  hard 
sphere  cutoffs  of  the  potentials  are  not  additive. 

6.  Conclusions 

A  combined  neutron  and  X  ray  diffraction  study  has  been  per 
formed  for  melt  spun  glasses  from  the  Ca6oMg4o-xCux  bulk  metallic 
glass  forming  system.  The  correlation  functions  show  clear  evidence 
for  solute— solute  nearest  neighbour  interactions,  unlike  some 
diffraction  studies  of  ternary  bulk  metallic  glasses  in  the  literature 
with  lower  total  solute  concentrations.  A  simple  approach  is  devel 
oped  to  estimate  the  maximum  solute  atom  fraction  where  sol 
ute— solute  contact  becomes  topologically  necessary.  The  Cu— Mg 
and  Cu— Ca  distances  are  consistent  with  the  sum  of  covalent  radii, 
whereas  all  other  interatomic  distances  are  consistent  with  the  sum 
of  metallic  radii.  Simultaneous  fitting  of  multiple  peaks  to  both 
neutron  and  X  ray  correlation  functions  can  be  used  to  gather 
limited  amounts  of  information  on  the  structure  of  ternary  bulk 
metallic  glasses.  Only  the  coordination  numbers  determined  for 
Cu— Cu  and  Cu— Mg  correlations  are  of  reasonable  reliability,  due  to 
peak  overlaps  for  the  other  correlations  which  occur  at  longer 
distances.  The  Percus— Yevick  (PY)  approximation  for  binary  hard 
sphere  systems  has  been  used  to  critically  examine  the  results.  This 
shows  that  the  fit  results  for  copper  are  a  reasonable  description  for 
atoms  which  are  in  close  contact  but,  the  simulations  indicate  that 
the  distributions  of  interatomic  distances  include  a  broader  trailing 
edge,  due  to  atoms  which  are  almost  touching,  which  is  not 
described  by  the  fitting  of  a  symmetric  peak.  A  combination  of  peak 
fitting  and  PY  has  demonstrated  that  copper  has  ~  6.5  neighbours 
which  are  in  close  contact,  but  the  PY  simulation  for  Ca6oCu40  indi 
cates  that  the  total  copper  coordination  number  is  about  12—13. 
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